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Abstract 

Dibaryon candidates with strangeness S = —2, —3, —4, —5, —6 are studied in terms of the ex- 
tended quark delocahzation and color screening model. The results show that there are only a few 
promising low lying dibaryon states: The H and di-0 may be marginally strong interaction stable 
but model uncertainties are too large to allow any definitive statement. The SIJ = —3,1/2,2 
NO state is 62MeV lower than the NJl threshold and 24MeV lower than the AStt threshold. It 
might appear as a narrow dibaryon resonance and be detectable in the RHIC detector through 
the reconstruction of the vertex mass of the AS two body decay. The effects of explicit K and 
rj meson exchange have been studied and found to be negligible in this model. The mechanisms 
of effective intermediate range attraction, a meson exchange and kinetic energy reduction due to 
quark delocahzation are discussed. 

PACS numbers: 12.39.-x, 14.20.Pt, 13.75.Cs 
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I. INTRODUCTION 



One of the most remarkable achievements of theoretical physics in the past thirty years 
is the establishment and development of the fundamental theory of the strong interaction 
- quantum chromodynamics (QCD). Perturbative QCD has been verified by high energy 
experiments. However, the low energy physics of QCD, such as hadron structure, hadron 
interactions and the structure of exotic quark-gluon systems, are much harder to calculate 
directly from QCD. One needs effective theories and phenomeno logical models in these cases. 

Borrowing the idea of quasiparticles from condensed matter and nuclear physics, one can 
approximately transform the complicated interactions between current quarks into dynamic 
properties of quasiparticles and what is left to be studied are the residual interactions be- 
tween quasiparticles. One of the quasiparticles in QCD is the constituent quark. How to 
dress the current quark to be a constituent quark still poses a theoretical challenge; various 
effective theories Uj, have been developed to derive constituent quarks from QCD. The 
common point of view is that the dynamical generation of the constituent quark mass is 
closely related to spontaneous chiral symmetry breaking initiated by the formation of a qq 
condensate in the QCD vacuum. 

The constituent quark model has been quite successful in understanding hadron spec- 
troscopy and hadron interactions even though we have not yet derived the constituent quark 
model directly from QCD. There are also various versions of the constituent quark model 
based on different effective degrees of freedom. De Rujula, Georgi and Glashowj:4=| first put 
forward a quark-gluon coupling model based on constituent quark and gluon effective de- 
grees of freedom. Isgur obtained a good description of hadron spectroscopy based on this 

n 

model|5]. However, extension of the model to baryon interactions does not reproduce the 
nucleon-nucleon intermediate and long range interaction. 

One modification studied is the addition of Goldstone boson exchange on the quark 
levelBBBa, tKe sKo«.a.ge pa« of the .te.a.io. . desCbed by tKe ,ua..- 

gluon degrees of freedom and the medium and long-range parts are attributed to meson- 
exchange. This quark gluon-meson exchange hybrid model achieves a quantitative fit of the 
nucleon-nucleon (NN) and nucleon-hyperon (NY) scattering data. 

A different modification of the De Rujula-Georgi-Glashow-Isgur (RGGI) model, the quark 
delocalization and color screening model (QDCSM)|l^, has also been suggested. It main- 
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tains the Isgur Hamiltonian for single hadrons but modifies it for baryon-baryon (BB) in- 
teractions with two new ingredients. 

First, the two center single quark orbital wave function (WF) used in the quark cluster 
model is replaced by a delocalized quark orbital WF. The introduction of quark delocaliza- 
tion can be viewed either as taking into account the contribution of excited configurations 
or the distortion of each individual baryon due to their mutual interaction. This straight- 
forward method enlarges the variational Hilbert space. Its advantage is that it permits the 
six-quark system to choose a more favorable configuration through its own dynamics, while 
maintaining a tolerable level of computational complexity. 

Second, a different parametrization of the confinement interaction is assumed, in which 
the usual quadratic confinement and a color-screened quadratic confinement are used to 
parameterize the two body matrix elements of different quark orbits. The introduction of 
this parametrization is aimed at taking into account some of the nonlinear, nonperturbative 
properties of QCD which can not be described by two-body quark-quark interactions in 
multiquark systems, such as the formation of color fiux tubes connecting many quarks, the 
three gluon interaction and the three body instanton interaction. 

The QDCSM not only justifies the view of a nucleus as, approximately, a collection of 
nucleons rather than a single big bag with 3A quarks, but it also explains the long standing 
fact that the nuclear force and molecular forces are similar except for the obvious energy and 
length scale differences. It is also the model which re quir es the fewest adjustable parameters 
to fit the existing BB interaction datafl ll.QQQHQ. 

There have been debates on which constituent quark model and which effective degrees of 
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2l| . In our 



freedom are best to use for hadron structure and interaction studies [ 
phenomenological study of BB interactions with three different constituent quark models, we 
found that even though the QDCSM and the other two models appear to be quite different, 
they give similar BB interactions in 44 of the 64 lowest BB channels consisting of octet and 
decuplet baryons IJ]. A preliminary analysis of the origin of this surprising similarity has 
been produced[22|. This result also implies that quark delocalization and color screening, 
working together, do provide the intermediate range attraction described by meson exchange 
in other models. On the other hand, the different models do give characteristically different 
results in some channels. For example, the QDCSM predicts a strong attraction in the 1=0, 
J=3 channel, producing a dibaryon resonance, d*, in this channel, while the quark gluon- 
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meson exchange hybrid model predicts a strong attraction in the 1=0, J=0 channel, 
implying the existence of a strong interaction stable di-fl. One cannot expect scattering 
data to become available in these channels to test these model predictions, but the dibaryon 
states should be detectable to provide a check on whether these model predictions are 
realistic. 

The study of dibaryon states not only checks constituent quark models but also searches 
for new hadronic matter. The i7-particle has been assumed to be a six quark system from 
the very beginningj2^ and has been both a theoretical and experimental topic for a long 
time. The S=0, = 0~ d' dibaryon was assumed to be an NNvr s yste m and was a hot topic 
in the 1990's[liQ- We showed that the 5 = 0, / = 0, J = 3 d*^, 14, 15, 26, 27, 3 
is a tightly bound six quark system rather than a loosely bound nuclear-like system of two 
A's. An S = —3,1 = 1/2, J = 2 NQ state was proposed as a high strangeness dibaryon 
candidate {29! . Kopeliovich predicted high strangeness dibaryons, such as the di-Q with 
S = —6, using the flavor SU(3) Skyrmion model j^O], and Zhang suggested searching for the 
di-Q in ultrarelativistic heavy ion collisions^]. Lomon predicted a deuteron-like dibaryon 
resonance using R-matrix theory [^3] and measurements at Saclay seem to offer experimental 
support for its existence. 

In this paper, we calculate promising dibaryons with strangeness S = —2, —3, —4, —5 
using the extended quark-delocalization color-screening model to provide another reference 
spectrum for strange dibaryons and a possible further check of constituent quark models. 
The SIJ=003 d* and -600 di-fl dibaryons, which we considered previously, are included in the 
discussion to demonstrate the differences between the predictions of the extended QDCSM 
and other quark models. 

The extended QDCSM is briefly introduced in Section II. In Section III, we present our 
results. We discuss these results further in Sections IV and V, and conclude in Section VI. 
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II. BRIEF DESCRIPTION OF THE EXTENDED QDCSM 

The QDCSM was put forward in the early 90's. Details can be found in Refs. jiol 
. Although the intermediate range attraction of the A^A^ interaction is reproduced by 
the combination of quark delocalization and color screening, the effect of the long-range 
pion tail is missing in the QDCSM. Recently, the extended QDCSM was developed jl^, 
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which incorporates this long-range tail by adding vr-exchange with a short-range cutoff. The 
extended QDCSM not only reproduces the properties of the deuteron well, but also improves 
agreement with NN scattering data as compared to previous workjiol. 

The Hamiltonian of the extended QDCSM, wave functions and the necessary equations 
used in the current calculation are given below. Here we do not take into account the effect of 
any tensor forces. The details of the resonating-group method (RGM) have been presented 
in Refs.Q, 3- 

The Hamiltonian for the 3-quark system is the same as the well known quark potential 
model, the Isgur model. For the six-quark system, we assume 

.2 6 
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-a^Xi ■ Xj 



if i,j occur in the same baryon orbit, 
if i,j occur in different baryon orbits. 



rij < ro, 
otherwise, 



where ro is the short range cutoff for pion exchange between quarks. All the symbols 
have their usual meanings, and the confinement potential Vconf{rij) has been discussed in 

Refs.0Q- 

The pion potential, V^(rjj), affects only the u and d quarks. We take these to have a 
common mass, rriq = rrid = rriu, ignoring isospin breaking effects as they are small on the 
scale of interest here. 

The quark wave function in a given nucleon (orbit) relative to a reference center (defined 
by S) is taken to have a Gaussian form characterized by a size parameter, b, 



(f - S) 



3/4 



(2) 



The values of rUg, rris, b, and Oc are determined by reproducing the A — N mass 
difference, the nucleon mass, a hyperon mass and by requiring a stability condition. The 
quark-pion coupling constant gggT, is obtained from the nucleon-pion coupling constant by a 



slight (< 10%) correction to the classic symmetry relation, viz., 

^ = (M^/m,)^g)'|ie-^^V2 (3) 

where M^v is the nucleon mass and the last factor provides the correction due to the extent 
of the quark wavefunction in the nucleon. The color screening parameter, /i, has been 
determined by matching our calculation to the mass of the deuteron. All of the parameters 
are listed in Table I. 

Table I: Model Parameters 



mg,ms{MeV) b{fm) ac{MeV ■ fm ^) 


£8 
47r 


ro(/m) /i(/m 2) 


313,634 0.6015 25.14 1.5585 


0.5926 


0.8 0.85 


The model masses of all octet and decuplet baryons 


are listed 


in Table II. 


Table II: Single Baryon Masses in Units of MeV 






N S A S A 






theor. 939.0 1210.6 1113.6 1350.0 1232.0 
expt. 939 1193 1116 1318 1232 


1358.1 
1385 


1497.5 1650.1 
1533 1672 



We use the resonating group method to carry out a dynamical calculation. Introducing 
Gaussian functions with different reference centers Si i=l...n, which play the role of gener- 
ating coordinates in this formalism, to expand the relative motion wave function of the two 
quark clusters, we have 
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In principle, any set of base wave functions can be used to expand the relative motion 
wave function. The choice of a Gaussian with the same size parameter, b, as the single 
quark wave function given in Eq.Q, however, allows us to rewrite the resonating group 
wave function as a product of single quark wave functions; (see Eq.()H) below). This cluster 
wave function (physical basis) can be expressed in terms of the symmetry basis, classified 
by the symmetry properties, in a group chain which in turn allows the use of group th eory 
methods to simplify the calculation of the matrix elements of the six quark Hamiltonian 35| . 
In our calculations, we typically use 12 Gaussian functions to expand the relative motion 
wave function over the range 0-8 fm. For a few channels, such as the deuteron and H 
particle, 20 Gaussian functions are used to extend the range to 12 fm. 



After including the wave function for the center-of-mass motion, the ansatz for the two- 
cluster wave function used in the RGM can be written as 

3 6 



k i=l a=l p=A 



where k is the channel index. For example, for SIJ = —2, 0, 0, we have k = 1,2, 3, corre- 
sponding to the channels AA, NS and SS. 

The delocalized orbital wavefunctions, ipR^r, Si, e) and '^/'^(r, Si,e), are given by 

Mr, 5., = ]^ Ur + y ) + e<Pir- y ) 1 , (5) 



iV(e) = VI + + 2ee-^?/^*^ 



where ^) and are the single-particle Gaussian quark wave functions referred to 

above in Eq.(0), with different reference centers y and — ^, respectively. The delocalization 
parameter, e, is determined by the dynamics of the quark system rather than being treated 
as an adjustable parameter. 
The initial RGM equation is 

J H{R, R')x{R')dR' = E J N{R, R')x{R')dR'- (6) 

With the above ansatz, the RGM Eq.® is converted into an algebraic eigenvalue equation, 

ECj,k'H^f = EY.C,,kNl^, (7) 

j,k' j 

where N^-,H^'^ are the wave function overlaps and Hamiltonian matrix elements, respec- 
tively, obtained for the wave functions of Eq.(j3]). 



III. RESULTS AND DISCUSSION 

Previously, we chose the di-f^ as an example to study whether or not our model results 
were sensitive to the meson-exchange cut-off parameter, tq, and the result demonstrates that 
they are notjs^. Hence, we consider it sufficient to calculate six-quark systems of different 
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strangeness, 5* = —2, —3, —4, —5, with a representative cutoff value of tq = 0.8/m. Table 
III displays the masses (in MeV) calculated for the strange dibaryon states of interest here. 
The lowest (without taking tensor coupling into account) channel for each SIJ is identified 
by bold lettering; sc and cc denote single channel (the lowest one) and coupled channels, 
respectively. '*' denotes those states without coupled channels. Additionally, it should be 
noted that in our calculation we assume the wavefunction to be zero at the boundary point, 
which is the usual boundary condition for bound states. If the state is unbound, we will not 
obtain a stable minimum eigenenergy in the course of extending the boundary point. The 
unbound states are denoted by "-" in Table III. 

Table III : Masses of Six-Quark Systems with Strangeness 
S,I,J coupling channels Massac MasScc 



-2,0,0 AA-iVS-SS - 2225.5 

-3,1/2,2 Nn - SS* - SS* - S*A - 2566.4 2549.1 

-3,1/2, 1 AH-iVr]-AH*-SS*-S*S*-SH*-SS - 

-4,1,0 H3-S*fi-S*S* 

-4,0,1 EE - AQ - EE* - 

-5,1/2,0 E*n 3145.0 * 

-5,1/2,1 En-E*n 

-6,0,0 3298.2 * 



In 1977, Jaffe|23l] studied the color-magnetic interaction of the one-gluon exchange po- 
tential in the multiquark system and found that the most attractive channel is the flavor 
singlet with quark content u^(Ps'^ . Moreover, the same symmetry analysis of the chiral boson 
exchange potential also leads to the very same conclusion 

However dibaryon physics can be very delicate. The deuteron channel is not a channel 
with strong attraction in any baryon interaction model. If the deuteron had not been found 
experimentally, it seems highly unlikely that any model would have been able to predict it to 
be a stable dibaryon. The if -particle {SIJ = —200) is a six quark state consisting mainly of 
octet-baryons, similar to the deuteron, and we find only a weak attraction there in our model 
also. Hence, a qualitative analysis is insufficient to judge whether or not the if-particle is 
strong interaction stable. Systematically, we find that a strong attraction develops only in 
decuplet-decuplet channels and a mild attraction in octet-decuplet channels. 
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Moreover, in the i7-particle case, the channel couphng effect may even be more important 
than the deuteron case. In fact, it is not bound without taking coupled channels into account. 
In our calculation, three channels have been taken into account. These are: EE, NS and AA. 
The relative motion wave functions of each channel contribution to the /J-particle are shown 
in Fig 1. We find that the A A channel provides the largest contribution (67%), followed by 
the NS channel (23%); the SS channel contributes only (10%). However, according to the 
analysis by Jaffe, the biggest contribution is the NS channel, and the AA channel provides 
the smallest contribution. 

The lowest mass we find for the u^d^s^ system is 2225.5MeV, which is 6MeV lower 
than the experimental threshold of AA and 1.7MeV lower than our model threshold. These 
values are smaller than the rms uncertainty that may be inferred from our fit to the baryon 
octet and decuplet in Table II. Furthermore, these values are on the order of corrections one 
would expect from isospin violations which we have not included. Hence, we can draw no 
definite conclusion as to whether or not the if-particle is strong interaction stable in our 

n 

model and we consider this to be consistent with recent experimental findings j87|. 

Besides the binding energy of the H, an interesting question regarding the H is its 
compactness, i.e., whether the H is a compact 6-quark object or a loosely bound AA state. 
Fig.l indicates that the maxima of the relative motion wave function of the dominant AA 
channel occurs around 1.8 fm, and the delocalization parameter, e, of the dominant channel 
at the maximum is ~0.1; hence, the if is a loosely bound system similar to the deuteron in 
our model. Such a similarity may well be more physically sensible than a compact 6-quark 
structure. 

For systems with strangeness S = —3, we have calculated the state Nf2 (SIJ=-3,l/2,2), 
which was shown to be mildly attractive, with energy below ASvr threshold j29|. That con- 

and supported by Silvestre-Brac and Leandri 



elusion was challenged by Oka 

We have carried out a dynamical channel coupling calculation to examine this state 
further. The NQ, AS*, SS*, SS*, channels are all included. We find the eigenenergy to 
be 2549.1 MeV, 24(54) Me lower than the AHtt experimental(model) threshold. Relative 
motion cluster wave functions for the individual channels are shown in Fig. 2. The NQ 
channel is by far the dominant one (77%) and the maximum for its relative motion wave 
function occurs at around 0.8 fm. The value of e at this separation is 0.46. Mixing into the 
other channels is small. Hence, we find this to be a compact six quark state. 
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Note that the D-wave AS and ES channels have not been included since no tensor inter- 
action has been included in this calculation. This coupling should be weak because there is 
no Trexchange in these channels and K and rj exchanges are not important in our model (see 
next section). Its effect on the eigenenergy of the NQ (SIJ=-3,l/2,2) state should be small 
and the D-wave decay widths to AS and SS final states should be small also. Therefore 
we expect the SIJ — —3, 1/2,2 Nfl should be a narrow dibaryon resonance. The tensor 
coupling calculation is currently in progress. 

We have also calculated the state SIJ — —3,1/2,1. In the course of extending the 
boundary point, the lowest eigenenergy closely approaches the model threshold but fails 
to come to a stable value within our limits of computation for the extension. This result 
is unchanged by taking into account all possible coupling channels. Due to the fact that 
there are only weak attractions in our model for the octet-octet channels and the size of 
the model uncertainty, it is difficult to conclude whether or not there are strong interaction 
stable states in these channels. We can only conclude that we have not found evidence for 
a strong interaction stable state with SIJ — —3, 1/2, 1. 

For systems with S=-4, we take the quantum numbers SIJ = —4, 1, as an example, 
because this case bears a number of similarities to the dcutcron. For example, in both cases, 
the lowest mass channel is composed of two octet baryons from the same isodoublet. Also, 
the matrix element P^/'^ characterizing the symmetry property of the system is -1/81 for 
both cases. (P3/ is the permutation operator of the quarks between two clusters acting 
in spin, flavor and color space.) The result shows that the system with S=-4, 1=1, J=0 is 
unbound, even when the S*S* and S*Q channel couplings are taken into account. Because 
the calculated energy is again very close to the SS threshold this conclusion should be viewed 
as tentative. 

For comparison, we have also calculated the SIJ — —4, 0, 1 state as shown in Table III. 
The SS, SS*, AQ and S*S* couphng channels are included. Our model result is very similar 
to the SIJ — —4, 1, 0, i.e., we do not find a bound state in this channel. The ^Si Di 
tensor coupling should be small in the SS channel so taking into account the tensor coupling 
will not change the unbound character. 

For the system with S=-5, we take the SIJ = —5, 1/2, state as an example. This state is 
interesting due to its < P^q'^ > value, ~ —1/9, which makes it a Pauli principle favored state. 
If only two-baryon S-wave channels are taken into account, there is only one channel for this 
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state. Our calculation shows that the contribution of the kinetic energy term, due to quark 
exchange and delocalization effects, contributes strongly towards the formation of a bound 
state. However, the one-gluon-exchange interaction largely compensates for this attraction 
and produces a mass of 3145. OMeV, which is 59MeV lower than the experimental value of 
the threshold but only about 2MeV lower than the model threshold. (The one-gluon- 
exchange effect here is quite different from that in the d* case, where large delocalization 
is favored for a wide range of cluster separations as well as there being a strong effective 
attraction due to the large reduction in the kinetic energy that accompanies significant 
delocalization.) We conclude that this state is not a good candidate for a dibaryon resonance 
search due to its small binding and its content. 

Inclusion of the tensor interaction will mix the spin-2 D-wave channel with the spin-0 
S-wave S*n channel. The tensor coupling should be weak also and we expect that its effect 
on the SIJ = —5, 1/2, state is small. 

In the same strangeness sector, we also calculated the SIJ = —5,1/2,1 state, since it 
includes the lowest channel, SQ. The calculated energy is very close to the threshold 
but a little higher so again there is no bound state, with SIJ = —5, 1/2, 1, in our model. 

To sum up, there are only a few high strangeness states worthy of experimental searches 
in our model. These are the H particle, the NQ and the di-Q. The di-Q was previously 
reported in Ref.( 3^) and the result is included in Table III. 

The if-particle and the di-Q may be strong interaction stable. However in our model, 
the binding energies of both are small relative to the model uncertainty. The di-f2 mass is 
about A7MeV lower than the experimental Qfl threshold. However our model mass for the 
Q is 1650MeV. If this model mass of Q were used to calculate the threshold, then the di-Q 
mass is no more than 2MeV below that threshold. Since our model mass for the single Q 
baryon deviates from the experimental value about 22MeV^, a reasonable estimate of the 
model uncertainty for the dibaryon would be at least that large. Therefore the di-f2 should 
not be claimed as a strong interaction stable dibaryon within the model. Similarly, we 
cannot claim that the if-particle is strong interaction stable either. The SIJ = —3, 1/2,2 
NQ case is certainly not strong interaction stable. However, the state is also certainly lower 
in mass than the NQ threshold, and quite possibly lower than the AHvr threshold, as well. 
The tensor coupling to the AS and SS channels should be weak and the decay width should 
be small. This strongly suggests that it is a promising candidate for a narrow dibaryon 
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resonance. This prediction can be tested by relativistic heavy ion reactions using RHIC 
detectors through the reconstruction of the vertex mass of the two body decay products, A 
and S. 



IV. THE K,ri MESON EXCHANGE EFFECT IN THE EXTENDED QDCSM 

The effects of K and rj meson exchange have been studied in the di-f2 channel and found 
to be neghgible in our model j^. In this section, we carry out a further systematic study of 
the effect of {K, if) meson exchange on the masses of strange dibaryon candidates. A flavor- 
symmetric octet meson-quark coupling is assumed for all of the octet mesons (vr, rj). The 
quark-meson exchange potential has the usual pseudoscalar meson exchange form except for 
a short range cut-off, 

^ Arc 12m,m,- r,-,- 

a=4 I' J i-j 

47r IZminij Vij ■' 

A unified cutoff of tq = O.Sfm is used. This is aimed at avoiding double counting since, 
in our model approach^], short and intermediate range interactions have been accounted 
for by the combination of quark delocalization and color screening. The model parameters 
(6, ttc, as) are refitted as before, and are listed in Table IV. The recalculated masses of octet 
and decuplet baryons are hsted in Table V. 

Table IV: Model Parameters 



mg,ms{MeV) b{fm) 


a^iMeV ■ fm-^) 




ro{fm) 


/i(/m 2) 


313,634 0.6022 


25.03 1.5547 


0.5926 


0.8 


0.90 


Table V: Single Baryon M 


asses in Units of MeV 








N S 


AHA 


S* 




Q 


theor. 939.0 1217.47 


1116.90 1357.56 1232.0 


1359.61 


1499.70 


1652.27 


expt. 939 1193 


1116 1318 1232 


1385 


1533 


1672 



Comparing Tables I, II and IV, V, it is apparent that the addition of K and r] exchanges 
has modified the model parameters and single baryon masses only slightly. We also find that 
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the properties of the deuteron are reproduced as well as before by a minor readjustment of 



the color-screening parameter n 



36| 



With the full octet meson exchange, we recalculated every single channel case with dif- 
ferent quantum numbers. For the iZ-particle, we calculated both the single channel and 
coupled channels, because it is very sensitive to channel coupling. The results are presented 
in Tables VI and VII. For comparison, the results with vr exchange only are also listed. 

Table VI : Masses of Six-Quark Systems with n exchange only and full octet pseudoscalar 

meson exchange. The masses are given in MeV. 

S,I,J single channel only vr exchange vr, /T, exchange 

-3,1/2,2 Nn 2566.37 2556.95 

-3,1/2,1 AH - 

-4,1,0 EE - 

-4,0,1 EE - 

-5,1/2,1 EQ - 

-5,1/2,0 E*Q 3145.01 3146.15 

-6,0,0 nn 3298.20 3300.00 

Table VII : Mass of iJ-particle {SIJ = —2, 0, 0) with vr exchange only and full octet 
pseudoscalar meson exchange. The masses are given in MeV. 







only TT exchange 


7r,K,r] exchange 




EE 


2280.40 


2282.67 


single channels 


NE 


2263.95 


2268.12 




AA 






coupled channels EE - 


~NE- 


- AA 2225.48 


2230.28 



Tables VI and VII show that those systems which are unbound in the extended QDCSM 
with 71 exchange only, such as SIJ = —3, 1/2, 1; —4, 1, 0; —4, 0, 1 and —5, 1/2, 1, remain 
unbound after adding K and rj exchanges, while those systems, which are bound in Table 
III, remain bound. Moreover, their masses are almost unaffected by the addition of K and 
rj exchange. The largest difference is not more than lOMeV. Especially for the if-particle, 
the AA channel would be unbound without taking into account the NE and EE channel 
coupling; this character is not affected either by adding K and rj exchange. Even the AA, 
NE and EE channel mixing fractions, 56%, 23% and 21%, respectively, are very similar to 
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the case with vr exchange only. All of these results confirm our expectation that heavier 
meson [K, rj) exchange has already been mostly accounted for by the quark delocalization 
and color screening effects in our model approach. Hence, explicit inclusion of K and rj 
exchanges beyond the cutoff scale is not important in our approach. 



V. FURTHER DISCUSSION ABOUT THE MECHANISM OF INTERMEDIATE 
RANGE ATTRACTION 

We have reported that the QDCSM gives very similar effective BB interactions to other 
models ^ , in general. However a careful comparison found that our model results for the 
binding energy of high strangeness dibaryons are systematically smaller than those of the 
chiral quark model [40]. The difference is mainly due to the different mechanism for the 
effective BB intermediate range attraction. 



In the chiral quark modelj 



the intermediate range attraction is attributed to a 
meson exchange. Because of its scalar-isoscalar character, a meson exchange provides a 
universal attraction independent of the flavor. 

In the Bonn meson exchange model the a meson is an effective description of correlated 



his point has been confirmed by a chiral perturbation calculation of 



two vr exchange. 

the NN interaction 4^ . Such an effective a-N coupling can not be extended from the NN 
channel to other channels with strangeness by a universal cr-baryon coupling, because in the 
NN channel, there are NN, NA and AA intermediate states for two tt exchange, while for 
the AA channel, there is only the SS intermediate state, and for the QQ channel, no such 
intermediate state is possible. Therefore it is not justified to fix the parameters of a exchange 
in the NN channel and then directly extend that exchange to channels with strangeness. 

There are arguments based on spontaneous chiral symmetry breaking for introduction of 
a meson quark coupling. In the SU(2) case, the non-linear realization of chiral symmetry 
can be linearized and in turn the a and vr meson coupling constants can be unifiedj^, 0|. 
However one should note this universal coupling is restricted to the u and d quarks and 
the effective a is still due to even number multiple vr exchange. In the SU(3) case, one 
can introduce SU(3) chiral symmetry by neglecting the difference between the s and light 
quark masses, followed by spontaneous chiral symmetry breaking. However, the nonlinear 
realization of SU(3) can not be linearized in the same way as in the SU(2) case to obtain 
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a universal u, d and s quark a coupling. The internal quark structure of the a meson is 
a controversial issue, but an equal mix of uu + dd + ss is quite unexpected and so also is 
a universal u, d and s quark a coupling. Hence, a strong attraction in high strangeness 
channels arising from such a universal a quark coupling is quite questionable also. 

In the QDCSM, quark delocalization and color screening work together to provide ap- 
propriate short-range repulsion and intermediate-range attraction for different channels. We 
illustrate this mechanism by showing contributions of the kinetic energy, confinement, color 
Coulomb, and color magnetic terms to the effective BB potential Vbb{S), as well as the total 
sum, in curves a-e, respectively in Figs. 3-6 for a few typical channels. The value of e varies 
with the separation S and is also listed in Figs. 3-6. The contribution of tt,K,t] exchange 
with a cut-off tq = 0.8/m is small so we do not show it. These figures are ordered in terms 
of increasing values of < P^l'^ >, from Pauli favored to Pauli forbidden. 

The effective baryon-baryon interactions shown in these figures are obtained in the Born- 
Oppenheimer approximation. 



Vbb{S, 



{^bb{S,)\H\^bb{S,)) {^bb{S,)\H\^bb{S,)) 



i^BBiSdl'^BBiS,)) {^BBiS,)\^BBiS,)) 

Where the "^BB^Si) is the antisymmetric six quark cluster state at a specified separation. 
Si, as given in Eq. (4) without summation over k and i; H is the six quark Hamiltonian (1). 
The contribution of each term in Eq. (1) is defined similarly, for example. 



T/™«/^C^_ {'^BB{S^)\Vconf\^BB{S^)) (^BB (^i) | Kon/ 1 
^BB K^i) — 



{^bb{S^)\^bb{S^)) {'^bb{S^)\^bb{S,)) 
Quark exchange alone induces a weak reduction in kinetic energy. Quark delocalization 
enhances this kinetic energy reduction. The kinetic energy reduction is also dependent on 
the strangeness of the channels due to the inverse quark mass dependence of the quark 
kinetic energy. The higher the strangeness, the smaller the contribution of quark kinetic 
energy and hence the smaller the reduction of the kinetic energy due to delocalization (see 
Figs. 3-6 (a)). This makes the intermediate range attraction weaker for the higher strangeness 
channels and is the main reason that our model gives less binding for the S = —4, —5, —6 
states than does the chiral quark model|4o|. 

If the usual quadratic color confinement is used, it does not contribute to the effective 
BB interaction. With the introduction of quark delocalization, the usual color confinement 
contributes an effective repulsion, as shown in Fig. 7. The color Coulomb has a similar 
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behavior, as shown in Figs.3-6(c); it produces almost no contribution to the effective BB 
interaction without quark delocahzation, and contributes an effective repulsion with quark 
delocalization. These two terms working together almost totally forbid quark delocahzation. 
This implies that the internal structure of the baryon is unaffected by the mutual interaction. 
However, that is inconsistent with the observed difference between the nucleon structure 
function in a nucleus and in isolation as seen in deep inelastic lepton scattering (EMC 
effect). 

On the other hand, there is no compelling reason to assume that the two body confine- 
ment potential is a good approximation for a multi-quark system. At the very least, the 
three gluon and the three body instanton interactions, which do not contribute to the qq 
meson and q^ baryon but do contribute to the multi-quark system, have been omitted in 
every two body confinement model. The color flux scenario revealed in the lattice QCD 
calculation of two and three quark systems raises questions regarding an additive two body 
confinement approximation. To take these facts into account, the QDCSM reparametrizes 

j,yyyi 



the confinement by introducing color screening|lfll. Illl. ll2L ll5L l2fiL l28j . Figs.3-6(b) show that 
after the introduction of color screening, the confinement term contributes an additional 
attraction, which mainly reduces the repulsive core of the effective BB interactions (See 
Fig. 7). There is a systematic uncertainty in our model related to this term which is not yet 
quantified. 

The color magnetic term generally contributes to a repulsive core except for a few Pauli 
favored channels where it contributes an additional intermediate range attraction so that 
these channels develop a strong effective attraction. It is also dependent on the strangeness 
of the channels: The higher the channel strangeness, the weaker the color magnetic con- 
tribution, due to the inverse strange quark mass in the color magnetic term. Figs. 3-6 (d) 
display these results. 

Altogether, these figures show that in the QDCSM the intermediate range attraction is 
mainly due to reduction of the quark kinetic energy, and the degree of reduction of the quark 
kinetic energy is connected with the degree of quark delocalization; the latter is determined 
by the competition between these four terms and that competition is different in different 
channels. We have mentioned before that in the d* channel this competition produces a 
large quark delocalization for a wide range of separations between two A-like quark clusters 
so that the kinetic energy acquires a correspondingly large reduction, which in turn gives 
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rise to a very strong attraction in this channel. (See Fig.3). The SIJ = —6,0,0 has the 
same < P^^^ > as that of the SIJ = 0, 0, 3 channel and has a similar opportunity to develop 
large quark delocalizations and strong attraction. However, there the competition does not 
allow as large a quark delocalization to develop and so the resulting attraction is not as 
strong as in the d* channel. (See Fig.4.) 



VI. CONCLUSION 



To sum up, we have carried out a dynamical calculation for the most promising 
dibaryon candidates with high strangeness by using the extended QDCSM. Only a few 
high strangeness states recommend themselves for experimental searches. These are the H 
particle, the NQ and the di-Q. We recommend searching for the {SIJ = —3, 1/2,2) using 
NQ — * AS D-wave decay vertex mass reconstruction with the RHIC detectors as the case 
with the best likelihood. 

In the QDCSM, it is the quark dynamics that controls the competition among the four 
terms in the Hamiltonian (1): the kinetic energy, the confinement, the color Coulomb, and 
the color magnetic, and which determines the overall effective BB interactions. These are 



quite similar, though not identical, to those of the quark-gluon- meson hybrid model Il4| for 
the majority of BB channels. The exchange of explicit K and r] mesons is not important in 
our model for mass estimates of strange dibaryons because of our short range cut-off. 

This model, which has the fewest parameters, describes the properties of the deuteron 
and the existing NN, NA and NS scattering data. Up to now it is the only model which gives 
an explanation of the long standing fact that the nuclear and molecular forces are similar 
in character despite the obvious length and energy scale differences and that nuclei are well 
described as a collection of A nucleons rather than 3A quarks. In view of the fact that the 
//-particle has not been observed experimentally, the BB interaction in the AA channel 
predicted by this model may be a good approximation of the real world. Further refinement 
is possible by including more channel couplings and spin-orbit and tensor interactions. 

Of course, the QDCSM is only a model of QCD. The high strangeness dibaryon resonances 
may be a good venue for determining whether the QDCSM mechanism for the intermediate 
range attraction is more realistic than that of a universal u, d, s quark-cr meson coupling. 
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FIGURE CAPTIONS 

Fig.l Relative motion wave functions of the coupled channels for the SIJ=-2,0,0 state 
with eigenenergy 2225.5 MeV. 

Fig.2 Relative motion wavefunctions of the coupled channels for the SlJ=-3, 1/2,2 state 
with eigenenergy 2549.1 MeV. 

Fig. 3 Contributions of kinetic energy, confinement, color Coulomb and color magnetic 
terms (a-d respectively) to the effective potential and the total (e) for the SIJ=0,0,3 AA 
channel with P^q'^ — —1/9. In each subfigure, the dotted curve is for delocalization param- 
eter e=0.0, dashed curve for e=0.5, dashed-dotted curve for e — 1.0; the solid curve is for 
the self-consistent value of e determined by the system dynamics. 

Fig.4 The same as Fig.3 for the state SIJ = -6, 0, nn with < Pg/^ >= -1/9. 

Fig.5 The same as Fig.3 for the state SIJ = 0, 0, 1 NN with < P^i" >= -1/81. 

Fig.6 The same as Fig.3 for the state SIJ = -2, 0, AA with < P^i^ >= 0. 

Fig. 7 The contribution of conventional quadratic color confinement to the effective BB 
interaction. The sohd curve is for e=0.0, dotted curve for e—0.5 and the dashed curve for 
e=1.0. 
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FIG. 1: Relative motion wave functions of the coupled channels for the SIJ=-2,0,0 state with 
eigenenergy 2225.5 MeV. 
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FIG. 2: Relative motion wavefunctions of the coupled channels for the SIJ=-3, 1/2,2 state with 
eigenenergy 2549.1 MeV. 
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FIG. 3: Contributions of kinetic energy, confinement, color Coulomb and color magnetic terms 
(a-d respectively) to the effective potential and the total (e) for the SIJ=0,0,3 AA channel with 
Pgg*^ = —1/9. In each subfigurc, the dotted curve is for delocalization parameter e=0.0, dashed 
curve for e=0.5, dashed-dotted curve for e = 1.0; the solid curve is for the self-consistent value of 
e determined by the system dynamics. 
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FIG. 4: The same as Fig.3 for the state SIJ = -6, 0, 00 with < P^" > 




FIG. 5: The same as Fig.3 for the state SIJ = 0, 0, 1 NN with < P^f >= 
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FIG. 6: The same as Fig.3 for the state SIJ = -2, 0, AA with < Pgg >= 0. 
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FIG. 7: The contribution of the usual quadratic color confinement to the effective BB interaction 
with and without quark delocalization and color-screening. 
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